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Abstract. We consider random lozenge tilings on the triangular lattice with holes Qi, . . . , Qn in some fixed 
position. For each unit triangle not in a hole, consider the average orientation of the lozenge covering it. We 
show that the scaling limit of this discrete field is the electrostatic field obtained when regarding each hole Qi as 
an electrical charge of magnitude equal to the difference between the number of unit triangles of the two different 
orientations inside Qi. This is then restated in terms of random surfaces, yielding the result that the average 
over surfaces with prescribed height at the union of the boundaries of the holes is, in the scaling limit, a sum of 
helicoids. 



Introduction 

The study of correlations of holes in random tilings was launched by Fisher and Stephenson [10], who 
considered in particular the monomer- monomer correlation and obtained exact data suggesting rotational 
invariance in the scaling limit. Motivated by this, we studied correlations of a finite number of holes 
of various sizes and shapes in [1], [2] and [3], and found that in the scaling limit they are given by a 
multiplicative version of the superposition principle for energy in electrostatics. In this paper we consider 
the discrete field of average tile orientations in a random tiling with holes, and provide the proof of the 
fact that the electric field comes about as the scaling limit of this discrete field. This result was announced 
in [4]. 

Consider the unit equilateral triangular lattice drawn in the plane so that some of the lattice lines are 
vertical. The union of any two unit triangles that share an edge is called a lozenge. Inspired by Feynman's 
description of the reflection of light in [9, Ch. 2] and given the connections between lozenge tilings with 
holes and electrostatics discussed in [2] we construct a discrete vector field as follows. Let Ai,...,A„ 
be holes in the lattice whose boundaries are lattice triangles of even side-lengths. As one ranges over the 
set of lozenge tilings of the plane with these holes (a portion of such a tiling is illustrated in Figure 5.4(a)), 
any given left-pointing unit triangle e in the complement of the holes is covered by a lozenge having one 
of three possible orientations: pointing in the polar direction 0, the polar direction 27r/3, or the polar 
direction Att/3. Define F(e) to be the average of these orientations over all tilings (the precise definitions 
are given in the next section). 

Now apply a homothety of modulus to our lattice with holes, and "drag outward" the images of 
the holes and the image of e so that they shrink to points zi , . . . , z„ and respectively zq in the complex 
plane, as i? — > oo. 

Our main result can be phrased as follows. 
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where rjo is i/ie wnii vector ■pointing jrom Zi in the direction of zq and c\v{Q) denotes the difference between 
the number of right- and left-pointing unit triangles in the lattice region Q. 

The figure above shows on the left an instance of the field F in the case of two oppositely oriented holes 
of side two, and on the right the corresponding Coulomb field that it approaches in the sc;aling limit. 

In Section 5 (see Theorem 5.1) we phrase this result in terms of random surfaces. Lozenge tilings of the 
plane with holes lift to certain discrete multi-sheeted surfaces. The above result implies that the average 
over these surfaces (an instance of which is illustrated in Figure 5.5) converges in the scaling limit to a sum 
of helicoids (for lifting surfaces as in Figure 5.5, this limit surface is of the type depicted in Figure 5.7(b); 
see [8] for an overview on a variety of settings in which helicoids arise). The key to proving this is that we 
obtain explicit formulas for the placement probabilities of lozenges in the scaling limit. 

Phrased this way, our result is similar in spirit to the result of Cohn, Elkies and Propp [5] on the 
local statistics for random domino tilings of the Aztec diamond, where an explicit surface is found as the 
scaling limit of the normalized height function. Another related instance of an explicit limiting surface, 
concerning the placement probability of lozenges in a random tiling of a hexagon, was studied by Cohn, 
Larsen and Propp in [7]. In [6], Cohn, Kenyon and Propp show that the normalized average height function 
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corresponding to domino tilings of regions that scale to an arbitrary simply connected region converges 
to the unique function that minimizes a certain surface tension integral. These authors also conjectured 
expressions for the local probability densities of domino configurations in large domino tilings. These 
conjectures were proved and generalized by Sheffield in [14]. Random surfaces corresponding to random 
dimer coverings of periodic bipartite planar lattices were studied extensively by Kenyon, Okounkov and 
Sheffield in [13]. 

One notable difference between our result and the above mentioned ones is the presence of holes. This 
causes the lifting surfaces to be multi-sheeted. The holes are also responsible for creating the field we are 
studying — without them the field would be zero. 

It may also be worth noting that for instance in [7] the boundary is the boundary of a lattice hexagon, 
and remains of fixed size throughout the scaling process, while in our case the boundary consists of the 
union of finitely many boundaries of triangular holes, each of them shrinking to a point in the scaling limit. 
Related to this is the fact that, unlike the results from the literature mentioned above, we consider the 
average of un-normalized lifting surfaces — the average of the normalized ones is identically zero. 

1. Definitions and statement of results 

Let Qj be a finite union of lattice triangular holes of side two with mutually disjoint interiors, for 

j = 1, . . . ,n (even when Qj is not connected, we still view it as a single generalized hole, often referred 
to as a multihole). The joint correlation ui{Qi, . . . , Qn) is defined in [3] (and recalled in Section 5 of the 
present paper) by means of limits of tori; the definition is readily extended to allowing some of the Qj 's to 
be lozenges. It follows from that definition that if L is any possible lozenge location, the probability that 
L is occupied by a lozenge in a random lozenge tiling with holes at Qi, . . . , Qn is given by the ratio 

a)(L,Qi,...,Qn) 
a)((5i, . . . ,Q„) 

(see Section 5 for the details). It follows thus that the vector F(e) = F(e; Qi, . . . , Qn) described in the 
Introduction is given by 

w{Li,Qi,...,Qn) U){L2,Ql,. . . ,Qn) u{L3,Qi,...,Qn) 

F(e) = . . „ ^ X ei H rjr: 77^62 H 777; 77^^3, (1-1) 

co(Qi,...,Q„) u}{Qi,...,Qn) w(Qi,...,Q„) 

where Li, L2 and L3 are the lozenge locations containing the left-pointing unit triangle e and pointing in 
the 0, 27r/3 and 47r/3 polar directions, respectively, and the ej's are unit vectors pointing in the directions 
of the long diagonals of the Lj 's. 

Note that we can specify the location of any left- or right-pointing unit triangle (for short, we will 
call them left- and right-monomers) by indicating the location of the midpoint of its vertical side. The 
midpoints of the vertical sides of the unit triangles in our lattice can naturally be coordinatized by pairs 
of integers using a 60° coordinate system with axes pointing in the polar directions ±7r/3. 

Let Ea,b be the east-pointing lattice triangle of side 2 whose central monomer has coordinates (a, 6). 
Let Wa.b be the similarly defined west-pointing lattice triangle. We regard both of them as holes. 

For any q € Q and any strictly increasing list of integers a = (oi, . .. ,as) for which qai € Z and the 
gaj's (equivalently, the goi's) are mutually disjoint, i = 1, . . . ,s, define the multiholes and 

by ' 

-^a ~ Eai,qai U . . . U Ea^^qa^ 

For a hole Q in the lattice, let Q{x,y) stand for its translation by the vector {x,y) in our coordinate 
system. We say that an integer divides a rational number if it divides the numerator of a lowest terms 
representation of it. 

We can now give the precise statement of our main result. 
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Theorem 1.1. Let Xq^^ , . . . , xi^^ , y^^^ , Vm^ , Zq^^ , • • • , and Wq^^ , • • • , w'if'^ be sequences of integers 
so that limfl:^oo /R = Xi, linifl^oo y\^^ /R = Vi, lim_R_^oo z^^^ /R = Zj and \im.R-^oo Wj^'' /R = wj for 
< i < m and 1 < j < n. Assume the {xi, yi) 's and (zj, Wj) 's are all distinct. 
Then for any multiholes E^^ , • • • , E^^ and W^^ , • • • , Wj^^ with 3| 1 — g, the field 



Fx, 



defined by (1.1) has orthogonal projections on our coordinate axes with asymptotics 

2(.xo - Xj) +yQ- yi 

' (xo - XiY + {xo - Xi){yQ - yi) + (yo - ViY 



( (R) {R.)\ 3 I 2(.xo - X,) +ya- ?a 



.1=1 



.\-t 2{xo-Zj)+yo-w, Uof-] (12) 

^ ^ (xo - Zj)^ + {xo - zj){yo - Wj) + {yo - wj)^ \ \R ' 



and 



p ( (R) {R)\ ^ ^ j V , XQ - Xj + 2{yo - y,) 

xo-Zj+2{yo-wi) 1+0^1^ (13) 

^^^ixo-zj)^ + {xo-zj){yo-wj) + {yo-wj)^j^ \rJ ' 

where s j and tj are the lengths of 34 and hj , respectively. 

Furthermore, for any e > and any bounded set B in the plane the implicit constants above are uniform 
over all choices of the limits for which each distance among the points {xo,yo), ■ ■ ■ ,{zn,Wn) G B is at 
least e. 

Note that in our oblique coordinate system the Euchdean distance between the points {x, y) and {x' , y') 
is i/(x — x')^ + (x — x'){y — y') + {y — y')'^, and the orthogonal projections of (x, y) on the coordinate axes 
arc x+^y and ^x + y. Note also that the multiholc E^,-^ 2 4 2s-2) consists of a contiguous horizontal string 
of right-pointing triangular holes of side two. Due to forced lozenges in its complement, E^^ ^ ^ 2s-2) 
precisely the same effect as the right-pointing triangular hole of side 2s that contains it; a similar statement 
holds for W-multiholes. Since ch{E^,) — 2si and ch(W^.) = —2tj, one sees that the theorem stated in the 
Introduction follows as a special case of Theorem 1.1. 

2. Reducing the problem to exact determinant evaluations 

One crucial ingredient for proving the results of [3] was an exact determinant formula for the joint 
correlation of an arbitrary collection of disjoint lattice-triangular holes of size two. The arguments presented 
there prove in fact a more general statement, which we will need in the current paper. 

Our determinant formula involves the coupling function P(x, y), x,y gZ specified by 

P{x,y) = ^f t-y-\-i-t)-^-'dt, x<-i (2.1) 



2TTi 



2-7ri/3 



and the symmetries P(x, y) = P{y, x) = P{—x — y — 1, x) (see Kenyon [12]), and the coefficients Ug of its 
asymptotic series 

00 

P(-3r - 1 -H a, -1 -h 6) ~ ^(3r)-'*-iJ7s(a, b), r ^ 00, a,beZ. (2.2) 

s=0 

Let r(a, b) and l{a, b) denote the right- and left-pointing monomers of coordinates (a, b), respectively. 



Proposition 2.1. Assum,e that {r{ai,bi), . . . ,r{am,bm),l{ci,di), . . . ,l{cn,dn)} can be partitioned into 
subsets of size two so that the monomers in each subset share at least one vertex. Then if m> n we have 

u){r{ai,bi), . . . ,r{am,bm),l{ci,di), . . . ,l{cn,dn)) = |det[Mp Mu]\, 



where 



and 



P{ai - ci,6i - di) 
P{a2 - ci,b2 - di) 

Mf 

P{am - ci,bm - di) 

Uo{auh + l) Uo{ai + l,bi) 
Uo{a2,b2 + l) Uo{a2 + l,b2) 

Uo{am,bm + i) Uo{am + i,bm) 



P{ai - Cn, bi - dn) 
P{a2 - Cn, b2 - dn) 



Um^_^{a-i_,bl + l) Um-n ^ (fll + 
t/m-n ^(09.69 + 1) Um-n ^(ao. + l,bo) 

+ 1) i7m-L_i(am + l,^m) 



Proof. Since the set of monomers can be partitioned into pairs, m + n is even; hence so is m — n. Thus 
we can prove the statement by induction on m — n as follows. To relate a given configuration to one for 
which m — n is 2 units smaller, bring in and extra hole, W(3r, 0), and let r — ^ oo. The details go through 
just as in the proof of [3, Proposition 3.2]. □ 

In order to obtain the asymptotics of the field F(e) we need to understand the three coefficients in (1.1) 
when the holes are Eg,^ , . . . , E^,^ , W\^^ , . . . , Wb„ ■ The asymptotics of the denominators in the special case 
x^-^^ = Rxi^ . . . ^w^^ = Rwn, Xi,. . . ^Wn G 3Z, is worked out in [3, Theorem 8.1]. Proposition 2.1 supplies 
exact determinant expressions for the numerators. Each is the determinant of a block matrix whose blocks 
consist in turn of (mostly) 2x2 blocks. To write them down it will be helpful to define the following five 
families of matrices. 

Given two lists a = (ai, . . . , as) and b = (6i, . . . , &t), let 



P{x—z-\-ai — bj — l,y—w-\-q{ai — bj) — l) 
P{x—z-\-ai — bj ,y—w-\-q{ai — bj) — 2) 

P{x—z-\-a2 — bj — l,y—w-\-q{a2 — bj) — l) 
P{x—z-\-a2 — bj ,y—w-\-q{a2 — bj) — 2) 



P{x—z-\-ai — bj — 2,y—w-\-q{ai — bj)) 
P{x—z-\-ai — bj — l,y—w-\-q{ai — bj) — l) 

P{x—z-\-a2 — bj — 2,y—w-\-q{a2 — bj)) 
P{x—z-\-a2 — bj — l,y—w-\-q{a2 — bj) — l) 



P{x— z-\-as—bj — l.y—w-\-q(as — bj) — l) 
P(x — z-\-as—bj ,y—w-\-q(as — bj) — 2) 



P{x— z-\-as—hj—2.y—w-\-q{as—bj)) 
P{x—z-\-as—bj'-l,y—w-\-q{as — bj) — l) 



(2.3) 



(the display shows the jth "column" of 2 x 2 blocks of Ax^y^z,w{^, t>); there is one for each j = 1, . . . ,t). 

5 



The second family of matrices is 



Ua(x+ai,y+qai) Uo(x+ai—l,y+qai+l) ■■■ Uk{x+ai ,y+qai) f/t (rr+ai — l.y+qai +1) 

Uo{x+ai + l,y+qai-l) Uo(x+ai,y+qai) ■■■ Ukix+ai + l,y+qai-l) Ukix+ai,y+qai) 

Ua(x+a2,y+qa2) UQ(x+a2 — l,y+qa2 + l) ■■■ Uk{x+a2,y+qa2) Uk{x+a2 — l,y+qa2 + l) 

Ug{x+a2 + l,y+qa2-l) Uoix+a2,y+qa2) ■■■ Ukix+a2 + l,y+qa2-l) Uk(x+a2,y+qa2) 



Uo{x+as,y+qas) Uo{x+as — l,y+qas + l) ■■■ Uk(x+as,y+qas) C/fc (aj+as — l,y+gaa +1) 
Uo{x+as+l,y+qas-l) Uo(x+as,y+qas) ■■■ Uk{x+as+l,y+qas-1-) Uk{x+ae,y+qas) 



The remaining famihes consist of one-row or one-column matrices. The third is 



(2.4) 



Cx,j/,^,ui(t>) — [ ••• P(x-z-bj,y-w-qbj-l) P(x-z-bj-l,y-w-qbj) ■■■] 



(2.5) 



The fourth family is 



^x,y,z,w{^^ — 



P{x — z + ai — l,y — w + qai) 
P{x - z + ai,y - w + qai - I) 

P{x - z + a2 - l,y - w + qa2) 
P{x - z + a2, y - w + qa2 - 1) 



P{x - z + as - l^y - w + qag) 
P{x - z + ag,y - w + qus - 1) 



(2.6) 



and the fifth 



Gx,y{k) = [Uo{x + l,y) Uo{x,y + l) ■■■ Uk{x + l,y) Uk{x,y + 1)]. 



(2.7) 



Let Li{x, y), L2{x, y), and i3(.T, y) be the the lozenge locations pointing in the polar directions 0, 27r/3, 
47r/3, respectively, and containing the left-monomer {x,y). Let Sj and tj be the lengths of the lists 
and hj, respectively, for i = 1, . . . , m, j = 1, . . . , n. Then Proposition 2.1 implies^ 



(xf\y[''^),...,W^^ (zi^Kwi^^)) 



detMi 



(2.8) 



^As in [3, Lemma 5.1], we use the fact that E{x,y) can be replaced by the union of the two right-monomers r{x — l,y) 
and r{x, y — 1) , and W(x, y) by ^(a: -|- 1, y) U l{x, y + I). 
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where 
Ml = 



P(0,0) C^(H)_^(H)_^<H)_^(K)(bl) ••• '^4H)_^(H)_^(^K)_^(H)(b„) G^(H) _^(«, (j/) 



-D (H) (fi) (R) (fl)(am) A (R) (R) (R) (ji)(am,bi) ••• A (r) (r) (r) (b) (a^, b„) -B (r) (r) (a^, J^) 

ii/m )'*'o 'i»0 '"^ )i/m j'^i i "^i '*'m iym i^n i "^n '*'m iym 

(2.9) 

and = ~ Z]j=i ~ 1- We set 5 = si H h and T = ti-\ h t„. 

Another application of Proposition 2.1 gives 

where^ 
Denote 

Throughout this paper ( = e^'^*/^. 

Proposition 2.2. If xf'^ /R ^ xq,..., w'i^^ /R w„ and Xq^^ = ao (mod 3), . . . , wif'^ = Sn (mod 3), we 

have 

1 1 detMf /1\ 

where 
M'l = 

n noio,l3o,li,5i(f\ ... nao,Po,ln,5n(f \ r;o:o,0o(O _ T' _ -1] 

f)ai,0i,ao,Po(„ \ Acti,Pi,'yi,5i(„ f\ ... ^ai,/3i,7n,i5n ('o, f \ R^i'/^ifsi .f — T — 1 "1 



(2.14) 



^ We write for the submatrix of A at the intersection of the rows with indices in I with the columns with indices in J; 
[n] stands for {!,..., n}. 
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and 



(2.15) 



wi; 



-l + (a-/3)-(7-5) Ca-l 



[z-x-{w-y)C,Y 



(j:;)(i-.CF-^ 



.l + (c-f3)-{'i-S) 



(j:;)(i-</cy-^ 



[z-x-(w-y)C,Y 



:-^+'°-^'-<^-^H.j-i)(i-gO' 

[z-x-(w-y)i]o+-^ 
[z-x-{w-y)C,]i+^ 



C-^+<°-'^'-'^-^'C+£7^)(i-9C)-+^-'' 

[z-x-{w-y)i\ =+o-i 
[z-x-(u)-j/)CJ=+3-i 



[z-x-(u;-y)f]J 



(j:;)(i-gc)^-^ 



-3+(c-/3)-(7-5)C J 



[z-x-(«;-y)C]^+i 
[z-x-(u;-y)f]J + i 



[z-a;-(tu-i/)Cl»+J-i 

^-l + (a-/3)-(T-a) ('»+j-2 



-(«'-y)CJ*+^- 



{there are t pairs of such columns, one for each j = 1,. . . ,t), 



(2.16) 



B^;yHs,k) = 



C-^+°-'^(g)(i-9C)° 



{x-yCY 



(x-J/C)" 



(x-ycr-^ 



(S)(i-9C)" 
(x-yC) '° 

c^+°-^(S)(i-gC)'' 

{x-y()~'' 



(x-yC)' 



{'•c-yC)~'° 

C-^+°-'^(g)(i-9C)° 
{x-yO~'' 



c^+°-"(.^)(i-gc)°-^\ /c-^+°-"(,°,)(i-gc)-A /c^+°-'^(,-i)(i-gcr-^\ /r^+^i^UOo-gcrii 



(x-yC)» 



(x-yC)« 



{x-yQ^ 



(x-yCy 



(2.17) 



.0+(..-g)-(T-^)(i_^^),--l 



-2 + (q-/3)-(7-5)|.j^_^^^3-1 

[z-x-(u)-y)C]J' 



(2.18) 



^ Here and throughout the rest of the paper (/(C)) stands for f{Q — /(C ^). 



x,y,z,w\'^ J 



-2+(c-/3)-(-,-5) 



[z-x-{w-y)C] 

^0 + (c-/3)-(t-«) (.■^_^^-|0 

|z-a;-(ui-y)Cl 



[Z-X-(U)-J/)C]2 

^0 + (c,-/3)-(7-5)(--^_^^-|l 

[2-a:-(u;-y)C]^ 



|,-2+(a 

[ 



-2:-(ii)-j/)Cl = 



anrf 



{x-ycy> 



\x-yCr I \{x-yO 



[z-x-{w-y)(]'' 



(2.19) 



(x-yC) 



{x-yC.)- 



(2.20) 



Furthermore, for any e > o?id ayiy bounded set B, the implicit constant in (2.13) is uniform for all choices 
of xo, ■ ■ ■ , Wn for which {xo,yo), . . . , (z„, w„) G B have all mutual distances at least e. 

Note that (2.13) (2.20) (plus 2 analogs) imply already, before evaluating the determinants, that in the 
scaling limit F(e) does not depend on the spacings between the side- two constituent holes of the multiholes 
E^. and W^. (although interestingly the independence of q is not yet apparent). 



Proof By (2.12), (2.8) and (2.10) we have 



Pi 



det Ml 



detM 



(2.21) 



It readily follows from the integral expression (2.1) for P that P(0,0) = 1/3. Regarding the first column 
of Ml as the sum of two vectors one of which is (1/3, 0, . . . , 0)-^ and using the linearity of the determinant 
yields 

1 det Ml 

- -I 
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Pi 



detM 



where 



In the special case x^^^ = Rx\, . . . , w'i^' = Rw„ and xq = ■ ■ ■ = Wn = (mod 3), Lemma 5.3 of [ec] 
shows how to write the main term in the asymptotics of det M as a determinant independent of R. This 
can be accomplished in the current case of general x[^^ , . . . , wi^^ by essentially the same arguments. 

Namely, consider the following operation on a square matrix X in which rows ii,. . . ,ik are of the form 
/(ci), . . . , f{ck), respectively, for some vector function /: Transform rows ii, . . . ,ik of X as 



Ml = 
„(«) 



-^l|(l,l)-entry set to 0- 



(2.22) 

(2.23) 



r/(ci)i 




/(C2) 




-f{Ck)- 





©°/(ci) 

(C2 - Cl)PV(ci) 



-(cfe - ci)(cfe - C2) . . . (cfc - Ck-i)V'' V(ci) 



(2.24) 



where V is Newton's divided difference operator, whose powers are defined inductively by 1?°/ = / and 

V-fic,) = (P-^- V(c,+i) - V-^fic,))/{cj+r - c,), r > 1. 
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This operation has an obvious analog for cohinins. 

As noted in [3, §5] (and as can be seen by looking at (2.11), (2.9), (2.3) and (2.4)), operation (2.24) 
can be applied a total of 2m + 2n different times to the matrix M: each row of block matrices in the 
expression for M given by (2.11) and (2.9) provides two opportunities (along the odd-indexed rows and 
along the even- indexed ones), and each column consisting of ^-blocks provides two more. Let M' be the 
matrix obtained from M after applying these 2m + 2n operations. 

Since operations (2.24) preserve the determinant (see [3, Lemma 5.2]), 

detM = detM'. (2.25) 

Moreover, by construction, the 2x2 blocks of M' are obtained by replacing Rxi ^ x[^\ . . . , Rwn ^ wi^'' 
in formulas [3, (5.6)-(5.9)]. Propositions 4.1 and 4.5 and the arguments in the proof of [3, Lemma 5.3] 
imply then that^ 

\2S m n 

^ n n («.fc-«.o'n n (^.^^-m' 

' j=l l<k<l<Sj j=l l<k<l<tj 



X det(M") iJ^'tSl<fc<i<m«'=«'+Sl<Jc<!<n*fc*'-II)r=lSr=l«'=*i}^ 



(2.26) 



where a*. = (a^i, . . . , afe^^ J and b; = (6a, . . . , hi^) for all k and I. 
By (2.25) and (2.26) we get 

2S m 



2-) n n ((^jk-ajifu n (^i^-M' 

^ j=l l<k<l<sj j=l l<k<l<tj 

X det(M") ij2{Ei<fc<i<,„si=s!+Ei<fc<i<„tfe*i-Er=iI^r=i«'=*i} 

+ O ^ij^^^l<'=<!<"> *'=*'+Sl<fc<i<n *fe*|-Slr=l ^1 = 1 _ (2.27) 



Propositions 4.1 and 4.5 imply that the implicit constant above is uniform in xi, . . . , w„. 

Next we turn to the asymptotics of detMi. Note that all the 2m + 2n operations of type (2.24) we 
applied to M are also well-defined as operations on Mi. Indeed, it is apparent from (2.9) and (2.3)-(2.6) 
that each of the 2m, + 2n times we applied (2.24), the extensions to Mi of the involved rows or columns of 
M are also of the form required by this operation. Let M[ be the matrix obtained from Mi after applying 
these 2m + 2n operations. 

A calculation similar to the one that gave (2.26) yields 

\ 2S-I-1 m n 

^) n n («.'^-%o'n n i^.k-b^ir 

' j=l l<k<l<Si j=l l<k<l<tj 



X det(M") _R^{^i<'!<i<m '*fc«!+IIi<fc<i<r, *'=*!-EfeLi E;r=i sfe*!}-!. (2.28) 
By the preservation of the determinant we get 

(-. \ 25-1-1 rn n 
^) n n i^ik-ajifxi n (v-m' 
^ i=l l<k<l<Sj j=l l<k<l<tj 

X det(Mf) ij^{^l<fe<i<m «fc«!+I3l<ie<!<n *'=*i-Er=l I3r=l 

_i_ Q ^fi^{i2i<k<i<m «'!*'+z]i<fc<i<„ **:*i-5Dr=i z;r=i sfctj-i^ ^ (2.29) 



*For a matrix A whose entries depend on a large parameter, [[AJJ stands for the matrix obtained from it by replacing 
each entry by the dominant part in its asymptotics as the parameter approaches infinity. 
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and the implicit constant is again uniform in .ti, . . . , w„ by Propositions 4.1 and 4.5. The statement of the 
Proposition follows now by (2.22), (2.27) and (2.29). □ 

The same analysis proves also the following result. 

Proposition 2.3. // x\^^ /R ^ xq,..., w'if^ /R Wn and x^^^ = ao (mod 3), . . . , u;i^^ = (5„ (mod 3), we 
have 



P2 ■■= 



1 1 detM^' /I 
'3 2TriR detM" '^"[r 



(2.30) 



where 



(7ao,/3o,7iii5i (f. '] 



f)ai,0i,ao,l3o(„ \ /ia!i,/3i,7i,i5i (-„ 4. \ 



^Xo,VO,Zn,W„ \^ri) 



G2t{S-T-l) 



iai,/3i,7n,i5n f \ B^i'/^ifsT S — T—}\ 



Aam,l3rr„Jn,Sn ( a f \ ftOirr^ Sm. I a Q — T — \\ 



^x,y,z,w\^) ~ 



and 



^-l + (<.-/3)-(T-a)(^_^g-)j-l \ / ^-3+(a,-/3)-(T-a)(^_^g-)j-l 

[z-x-(w-y)C]3 I \ {z-x-{w-y)i\i 



(2.31) 



(2.32) 



-l + c,-/3 



-3+a-/3 



(x-yC)-'' 



{x-yCr I \ {=o-yCr I \(x-K)-V \ (x-yQ-^ 

The implicit constant above has the same uniformity property as in Proposition 2.2. 



{x-vO~'' 



(2.33) 



3. Evaluation of the determinants and proof of Theorem 1.1 

The orthogonal projections of the vectors ei, e2 and of (1.1) on the x-axis of our oblique coordinate 
system are \/3/2, — \/3/2 and 0, respectively. Thus the orthogonal projection of F on the a;-axis is 



TP ^ 



(3.1) 



By (2.13) and (2.30), the part about of Theorem 1.1 will follow provided we show that (detMf 
AetMl^)/ detM" evaluates to the expression that together with (3.1), (2.13) and (2.30) imphes (1.2). 
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Proposition 3.1. The determinants of the matrices M'{ given by (2.14) and (2. 16) -(2.20), given by 
(2.31)-(2.33), and M" given by (2.15), satisfy 



detM{'-detM^' 
detM" 



Sk 



2{xa - Xk) + yo-yk 



.fe=i 



{xq - Xk)^ + {xq - Xk)iyo - yk) + (2/0 - Vkf 



2{xo - zi) +yo-wi 



^ (xo - zi)^ + {xo - zi){yo - wi) + (yo - wiY 



(3.2) 



A quick glance at (2.14) and (2.31) shows that M" and Mlj differ only in their first rows. Furthermore, 
by (2.18), (2.20), (2.32) and (2.33), the corresponding first row entries differ just by a factor of . This 
makes it very tempting to write the numerator on the left hand side of (3.2) as a single determinant, and 
try to use the method of factor exhaustion (after all, this was the method that proved successful in [ec] , 
where an explicit product expression for a special case of det M" is given; see [3, Theorem 8.1]). However, 
we were not able to prove (3.2) this way, despite getting frustratingly close (the only missing part was 
proving divisibility by one last type of linear polynomial divisor). 

The proof below derives (3.2) by a certain limiting process from a convenient specialization of [3, Theo- 
rem 8.1]. This has the advantage of being significantly shorter than a possible proof by factor exhaustion. 
Still, a solution of the latter kind would be interesting, as it should prove also a conjectured two-parameter 
generalization of (1.2) that works in particular with C, being an indeterminate. 

Proof. We first prove the statement under the simplifying assumption that ao = • • • = (5„ = 0. In 
this case Theorem 8.1 of [3] implies that for any indeterminates xi, . . . , Xm, yi, - ■ ■ ,ym, zi,. . . ,Zn and 
wi,. . . , Wn we have 



det M" = det 



/jo, 0,0,0 (q + ^ 

^xi,yi,zi,wi '■1/' 



iO,0,0,0 /„ f\ 



jo, 0,0,0 



.(Sl,in) 



iO,0,0,0 /„ f \ oO,0 



. (^rri) S ■ 



1). 



l<k<l<m 



l<k<l<n 



'[[Y[[{=o,-zi-ayk-wi)){x,-zi-^J^)y'''' 



k=l 1=1 



(3.3) 



Let e > 0. Note that the matrix in (3.3) depends on four lists of indeterminates — Xi, . . . ,Xm, J/i, • • • , J/m, 
Zi,...,Zn, and Wi,...,w„ — and two lists of positive integers si,...,Sm and fi,...,i„. Replace the 
indeterminates as follows: the x-hst by xo,Xi, . . . ,Xm, the y-list by yo,yi, ■ ■ ■ ,ym, the ^;-list by xq + 
e,zi, . . . ,Zn, and the w-list by yo,wi, . . . ,Wn- Replace also the s-list by 2, si, S2, . . . , Sm, and the t-list by 
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2,ti,t2, ■ ■ ■ ,tn- Then (3.3) provides an explicit product expression for the determinant of the matrix 



~) (^) <^xo%o'°i,'u'i(*l) 

~^ Cxo,yo,zi,wi{tl) 
Dxi,yi,xo+e,yo{si) Dxi,yi^xo+e,yo{Sl) ■^x'i,'yuZi,wi{^i'>ti) 



/^O, 0,0,0 U \ /=<0,0 (,.\ 

^xo,Vo,Zn,w„\''nJ '-^xom^ > 

Cxo,yo,Zn,Wn{^n) Gxo,yoi'^) 

■^x'uyuz„,w„ ' ^n) ^x'i,yi (*1 ' ^) 



.Dxm,ym,Xo + e,yoi^l) ^Xm,ym,Xo+e,yo{^l) ^'xmaim,Zl,Wli^mjtl) ■■■ ^Xm',Jfm,Z„,l«„ (•^'"' ,J/m ( ' '^i 

(3.4) 

where v = S -T -1, 



z-x-{w-y)C]^ / \ [z-x-{w-y)C]^ 



(3.5) 



lz-x-{w-y)C\ 
[z-x-{w-y)C] 

C"'(i-gC)\ 

[z-x-{w-y)Cl'-' 
[z-x-(w-y)C]'-' 



[z-x-(w~y)C]' 
[2-a:-(ui-j/)C] = 



(3.6) Dx,y,z,wis) = 



lz-x-{w-y)C\ 
[z-x-{w-y)C] 

C-'(i-gC)\ 

[z-x~{w-y)C\'-' 

C(l-qCV 
lz-x-{w-y)CY-' 



[z-x-(w-y)C]' 



[z-x-{w-y)C,]^ 



(3.7) 



and 



Gx,y{k) 



(x-yOO / \ (x-yO° / \ (x-yO 



iic-yC)- 



(x-yCj *'/ \ {x-yC) 



(3.8) 



The resulting product expression for detM^ has many common factors with the right hand side of (3.3), 
which simplify when taking detM^/ detM". After simplification we get 



detM, 



detM" 



\^{xo-xk-(iyo-yk))(x 

3 k=l 



'-^)Y'''[l[(^o+e-z,-C{yo-m)){^o+e-z,-?!Sl^)Y 



1=1 



,2 m 



fe=l 



(3.9) 
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Replace the first column in (3.4) by tlie negative of the sum. of the first two columns. Using = 1, 
-C~^ - = C~^) and -C - = 1) we obtain 



- det Me 



, I ^0,0.0,0 \ /^O, 0.0,0 /i \ /oO,0 

^ ~ i ~ ''^xo,yo,2i,'u>i{tl) ■■■ Cxo,yo,Zn,Wn{tn) Gxo,yo{^) 

AO, 0,0,0 (e \ F) ^ iO, 0,0,0 (q -f \ iO, 0,0,0 (a + \ rO,0 /„ ,,\ 

^xi,yi,xo+e,yo\^'^/ ^xi,yi,xo+e,yQy^l ) ■^xi,yi,zi,wi\^i' ''^J ''' ^xi,!/i,z„,u)„ V*! ' ■^xi,yi\^i': J 



A0,0,0,0 (q.) f) , Csi^l iO,0,0,0 C„ ... 40,0,0,0 /„ f \ dO,0 /„ 

(3.10) 

Denote the matrix on the right hand side above by A'' = (Afy)i<j,j<2S+i- As e — > 0, the only diverging 
entries of A'' are N\\ and A^22- It follows from (3.10) that 



detM - ^^"^ detiv''^+'i\^'''^ + ^ '"^ 



det7V[^^+'l\^'^ - detiv''^+'l\^'^ 

"^^-'^[2S+1]\{1} "*^^-'^[2S+l]\{2} 



+ 0(1), e^O, (3.11) 

~ ~ [25+ ll \ -f 2'1- 

where A'' is the matrix obtained from N by replacing Nn and A^22 by 0. But -^[2s+i]\{2} J^^* matrix 
M^' of (2.31), and A^fal+Jui^J^ = M" . Furthermore, by Lemma 3.2, det TVfll+Ji^^^^j^ = detMf . Therefore, 
(3.11) implies 

det Me 3 C"^ -C detMf -detM^' „ 

Extracting the coefficient of 1/e in the asymptotics of (3.9) as e — > and comparing it with the second 
term in (3.12) gives (3.2). 

Next we show how the case of general residues ao, . . . ,6n modulo 3 reduces to the above case of all zero 
residues. 

Consider the general matrix M" given by (2.14) and (2.16)-(2.20). Except for its first row and column, 
it consists of 2 x 2 blocks of the type in Lemma 3.3. Furthermore, the corresponding value of a is constant 
over each block ^"^'y^'^'J/'^' (sfe, ^/), and equals {a^ — 0k) — {li — Si). It is also constant over each block 
^xl:;yk {sk,S -T -1), and equals Uk - Pk- Thus in each 2x2 block of type (3.13) in M{' the value of the 
a of Lemma 3.3 is a difference of two quantities, the first being constant along each row of Mf , and the 
second constant along each column of M". Therefore Lemma 3.3 can be used to transform the general M[' 
matrix by row and column operations into the specialization of Mf when all residues are except ao and 
f3o, and the overall effect on the determinant is that it is multiplied by a, where a G {1,-1}. The very 

same row and column operations transform Mg' into Mj'Iq^^o s„=o and M" into M"|cj^=o,...,i5„=0) with 

the effect on their determinants being multiplication by the same a. Therefore 

detM('-detM^' _ det Mf |„i=o,...,5„=o - det Af;^'|aj=o,...,5„=o 



detM" detM"|„i=o 



.5„=0 



However, Lemma 3.4 provides determinant preserving operations that transform all three determinants 
above into their ao = 0, /3o = specializations. □ 
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Lemma 3.2. det ivj2f = detMf . 



Proof. Repeated application of Lemma 3.4 transforms the first row and column of the matrix -^[2s+i]'\{i} 
(see (3.5), (3.7) and (3.8)) into the first row and column of the matrix Mf (see (2.18)-(2.20)). All the other 
and det M" agree, and are left in agreement by the above applications of Lemma 3.4. 



entries of N, 



[2S+1]\{1} 
[2S+1]\{1} 



This proves the claim. 

Lemma 3.3. Let f be a function defined at ( and C~^; o.iT-d let a G Z. Let 



□ 



A{a) = 



(C+V(C)) (C-V(C)) 



Denote its rows by Ri and i?2, and its columns by C\ and C2. 

(a) . Simultaneously replacing {i?i ^ R2, R2 < -Ri — R2} turns matrix A{a) into A{a — 1). 

(b) . Simultaneously replacing {C2 ■*— Ci,Ci ■* Ci — C2} turns matrix A{a) into A{a + 1). 

Proof. Since C"^ = 1, the second row of A{a) is the same as the first row of A{a — 1). The negative of 
the sum of the entries in the first column of A{a) is 

- (r-V(c) - r"+V(r')) - (c+V(c) - c-'fiO) = cfio - rv(r') = icf iO) , 

as — C'^'*'^ = C^^^ for ^11 integers k. One similarly checks that the negative of the sum of the entries 
in the second column of A(a) equals (C°~^/(C))- This proves (a). Part (b) follows analogously. □ 
By a similar calculation one can easily check the following. 

Lemma 3.4. Let f be a function defined at ( and (^^, and let a, /3, 7 e Z. Form the matrix 







(C-'+"/(C)> 
(C-^+V(C)) 



Denote its rows by Ri, R2, R3, its columns by Ci,C2,C^. Then the result of the simultaneous column 

operations {C2 < C2— C3, C3 ^ C2}, followed by the simultaneous row operations {R2 < R2—R3, R3 <— 

R2}, is the matrix 

(C"/(C)) (C-'+"/(C))' 

<C-'+'^/(C)> (c-'^vio) (C-^+V(C)) 

. <c''/(c)> (e+vio) (c-^+v(c)). 



Proof of Theorem 1.1. Part (1.2) of the statement follows from (3.1), (2.13) and (3.2). Interchanging 
the roles of the coordinate axes we get (1.3). □ 



4. A finer integral asymptotics 

Let T> denote Newton's divided difference operator, whose powers are defined inductively by T>^ f = f 
and f{cj) — {V^~^ f{cj+i) — T>^~^f{cj))/{cj+r " Cj), r > 1. We will need the following result on the 
asymptotics of the coupling function P when acted on in the indicated way by powers of T>. 

Theorem 4.1. Let rn and s„ be integers so that lim„^oo fn/n = u, lim„^oo Sn/n = v, and {u, v) ^ (0, 0). 
Then for any integers k,l >0 and any rational number q with 3|1 — g we have 



T^'v {^l +X + y. Sr, + q{x + y)%=a, } 1^^,^ = 
2^i\ k 



{-rn 
15 



k+l+1 



o 



1 



(4.1) 



where C, = e^'^'^^^, acts with respect to a fixed integer sequence ai, a2, . . . , and Vy acts tuith respect to 
an integer sequence 61, 62, . . . satisfying qbj € Z for all j > 1- Furthermore, for any open set U containing 
the origin the implicit constant above is uniform for all {u, v) ^ U. 

We note that Proposition 7.1 of [3] corresponds to the special case when r„ = —un and s„ = —vn. It 
will be crucial in our proof of the scaling limit of the average lifting surface to know that the leading term 
on the right hand side of (4.1) is independent of the way r„/n and s„/n approach their limits. 

We deduce Theorem 4.1 from the following auxiliary results. Let V denote the counterclockwise oriented 
arc of the unit circle connecting ( = e^'^'/^ to —1; V[C, —1) stands for this path minus the point —1. 

Lemma 4.2. Let rn and Sn he integers, and assume lim„^oo rn/n = u> e> Q. Then if q € C^{V[C, —1)) 
has a pole of finite order at —1, for any integer k > we have 



(4.2) 



where the implicit constant depends only on k, e and q{t) {in particular, it works for all u > e and all 

integers .s„). 

Proof. Let / denote the integral above. Write q{t) = q{t)/{t + 1)', with / > and g G C°(P). Making 
the change of variable t = e'^ we obtain 



I = (-1)-" ^ \i + ty^-h'-it - o'^m 

= i-^Y- [ 
J 2 



dt 



2v/3 

tt/3 



2cos^e*^/2 



d9 







/3+r/2) 
k 



J(2-K/3+T)Sr, 



(4.3) 



Since the graph of 2cos(7r/3 + t/2) is concave for r G [0, 7r/3], it lies below its tangent at r = 0. This 
implies 

/IT T\ \/3 T 

'^°n3 + 2)^l-^^<l-2' 
The elementary inequality \e^ — 1| < \y\e^^^ implies 

|e*^ - 1| < re^ < 3r, rG[0,7r/3]. 

The above two inequalities combined with (4.3) give 



/|<_^^^'(l-0""'(3r)'=|g(e'(2-/3+^)) 



dr 



<3'=sup|g(i)| / (1--) 



r^dT. 



(4.4) 



However, integration by parts implies 

j-n/Z 

J{r,k) :-- 



T (-f)"' (5)"' 
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2k 



Repeated application of this, together with 

J(r + k,0) = r(l-^) + ; 

^ ^ ' ^ r+k+l\ 6/ r+k+2 

shows that J(r, k) is equal to a sum of fc + 1 terms each exponentially small in r, plus 

2 (2fe)(2fc~2)---2 
r + k + 2{r+ l)(r + 2) • • • (r + fc) ' 

Since r„ — / + 1, . . . ,r„ — Z + fc + 2 > ^nu for n large enough, it follows that the integral in the second 
line of (4.4) is majorized for n large enough by 2'^'^+^k\/{nu)''~^^ < 2^'=+^fc!/(e)*+^l/n''+^, and the proof is 
complete. □ 

Lemma 4.3. Let q{t) = {t-C,Yqi{t), where qi{t) G C\V), k is a non-negative integer, and gi(C) 7^ 0. Let 
Tn, Sn^'^ so that lim„^(x) Tnl'n = « > 0. 
(a). //fc = 0, 

/ (-1 - tY-t^-q{t) dt = +o(^). (4.5) 

{b). Ifk>l, 

f\-i - tyn^^qit) dt = ^r^^ - ty-H-m dt + o (^) . (4.6) 

Furthermore, for any e > 0, each implicit constant above can be chosen to be uniform for u e [e, 00) and 
independent of Sn- 

Proof For Z = 2, 3, . . . , 7, define 

^'W^=E7T^(*-C)'+'^ teV[C,-i). (4.7) 

i>0 •' 

One readily checks that the Taylor series expansions of In t and ln(— 1 — t) around t = ( can then be written 

as 

Int = InC + C^it - C) - C^Ih + C^h7 - Ch2 + (hn + /13 - he, 

and 

ln(-l -t)= ln(-l - C) - at - C) - (hi - Chr - C^h^ - - /13 - he- 

Therefore we have 

rln(-l-i) + slni= r ln(-l - C) + slnC + (sC"^ -<)(*- C) 

+ {-rC^ - s0h2it) + i-r + s)h3{t) + i-rC - sC^)h4{t) 

+ i-rC' + sQh^it) + i-r - s)he{t) + i-rC + sC-')hr{t), tGV[C-l). (4.8) 
Let r,s gZ and denote 

I{r,s) := j {-l-tYt'q{t)dt. 
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Using (4.8) and = 1, integration by parts gives 



/(r,s) 



rln(-l-t)+slnt 



q{t) dt 



-'-C)(t-<)g(-rC-i-sC)/»2+- + (-rf+sf-i)/i 



''q{t)dt 



«-i-rC)(t-C) 



where 



dt 



h{t) := {~rC^ - sQh^it) + (-r + s)/i3(t) + {-< " sr')/i4(i) 

+ {-rC^ + sQhr,{t) + i-r - s)he{t) + (-< + sC^)h7{t). 

Since C*~'^e(*f '^-'"C)(t-C)gb(t) _ (—1 — ^)''^^^ we see that the upper Hmit in the first term of the expression 
in the large curly braces above equals whenever r > k. Thus, for r > A; we obtain 



7(r,s) 



rC - sC 



M~'-rO{t-C)Mt) 



e^^*>q{t)h'2{t)dt + ' 



+ - 



-rC + 



-K)(t-f)g6(t) 



q{t)h'7{t) dt 



+ 



,K-i-rC)(t-f)^6(t) / 



e''^'>q'{t)dt. 



If Irl > Isl we have 



1 - s/rC 



< 



l + |s/r| 



< 



|Rc(l-,s/rC)| - 1/2 



4. 



(4.9) 



(4.10) 



A similar argument shows that the above inequality holds in fact also when |s| > \r\. All six fractions in 
front of the integrals in the expression in curly braces above are thus seen to be majorized in absolute value 
by 4. 

Regard /(r„, s„) as the sum of the three quantities provided by (4.9). To deduce part (a) of the Lemma, 
assume k = 0. All six terms of the second quantity are 0(l/n^) thanks to (4.10) and an application 
of Lemma 4.2 with k = 1 (which applies since h2{t), . . . ,h'^{t) arc all of the form (t — (^)g{t), where 
g e C°°{V[(, —1)) has a simple pole a.t t = —1), followed by an application of Lemma 3.4 with fc = for 
the resulting integral on the right hand side of (4.6). Finally, the third quantity provided by (4.9) is also 
0(l/n^), due to the fraction in front of the integral and another application of Lemma 4.2 with k = 0. 

For part (b), assume A; > 1. Then q{() = 0, and (4.9) provides an expression for I{rn,Sn) as a sum of 
just two quantities. All six terms in the first quantity are 0(l/n'^+^) due to (4.10) and Lemma 4.2 applied 
with k replaced by fc + 1 (as explained in the previous paragraph, this unit increment comes about by the 
presence of the h'i{t) factors in the integrands). This proves (4.6). 

The uniformity of the implicit constant follows because both the majorant in (4.10) and the implicit 
constant in (4.2) are uniform. □ 

Repeated application of part (b) of the above lemma and one final application of part (a) yields the 
following result. 
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Proposition 4.4. Let q{t) = (t - C)''9i(0; where k>0,qi€ C''+^(P[C, -1)) has a pole of finite order at 
t = —1, and gi(C) 7^ 0. Let rn, Sn so that lim„^oo rn/n = u> 0. Then 

Furthermore, for any e > the implicit constant is uniform for u G [e, 00) and independent of Sn- □ 

Proof of Theorem ^.L Suppose first that w < 0. Then (2.1) holds for large enough n, and for 3|1 — q 
[3, (7.11)] and [3, Lemma 7.4] give 

K {Dl Pirn +x + y,sn + q{x + y))U=a, } 1^^,^ 

k+l 



+ Ck+l+l{t-C)'+'+'+^ 



Proposition 4.4 applied to the right hand side above yields 

V'y {Vl P(r„ + X + y, s„ + q{x + y)%=a, } |^^,^ 

1 /C"'""fcfn(C-gr^)"+'(fc + 0!\ , 1 



which is just what (4.1) states. The uniformity of the implicit constant above follows by the uniformity of 
the implicit constant in (4.11). 

Since (u,w) 7^ (0,0), at least one of u < 0, f < 0, and —u — f < is true. The symmetries P{a,P) = 
P{—a — /3 — 1, a) and P{a, fi) = P{(3, a) of the coupling function allow one to use the same arguments that 
proved the case u < to deduce the other two cases (see the proof of Proposition 7.1 in [3] for details). □ 

Proposition 4.5. Let rn, Sn & so that lim„_>oo r„/n = u and lim„_>oo Sn/n = v. Then for any integers 
k,l>0 and any rational number q with 3\1 — q we have 

Vl Ui{rn +x,Sn + qx)U=a, = ^ Q (^''.-^.-1(1 _ qQ^^r^ - s^Ct') 

+ 0(n'-'=-i), (4.13) 

where acts with respect to some fixed integer sequence ai, 02, Given any bounded set B in the plane, 

the implicit constant can be chosen so that it is uniform for {u,v) G B. 

Proof By [3, (6.8)] one has 

Ui{ci,b) =l-{C-'-\a-bO') 

+ monomials in a and b of joint degree < I. 

This implies, for 3|1 — g, that 

Ui{rn+x,Sn + qx) = 7^(C''""'""M(l-9C)a; + (r„-s„C)]') 

+ Ca,l3{rn+xY{Sn+qxf, (4.14) 

a,/3>0 
a+/3<I 
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where Ca,p is independent of r„ and s„ for all a and /?. 

On the other hand, the argument that proved [3, Lemma 6.4] implies that for any constants A,B,CgC 

Vl{Ax + Brn + CsnYU^a, = A\Brn + CsJ-'^ + O {n'-'^-') , (4.15) 

with implicit constant uniform for {u,v) ranging over any bounded set. Combining (4.14) and (4.15) yields 
the statement of the proposition. □ 



5. Interpretation in terms of height functions 

Lozenge tilings of regions with no holes are well known to be interpretable as lattice surfaces (see e.g. 
[15] [7]). Regard the unit triangular lattice T on which the tiled region lives as being in a horizontal plane, 

and let £ be a copy of the lattice placed so that one family of its body-diagonals is vertical, 

intersecting this plane at the vertices of T. Then each segment joining two nearest neighbors of C projects 
onto a unit lattice segment of T. 

Orient the lattice segments of T so that they point in one of the polar directions 7r/2, — 7r/6, or — 57r/6. 
Then we can lift a lozenge tiling of a lattice region on T by starting from some lattice point, tracing around 
its tiles one after another, and at each traversal of a lattice segment ,s of T, moving cither up or down 
on the corresponding lattice segment of C, according as the traversal respected or violated the orientation 
of s. Tracing around a lozenge results in going around a lattice square of C whose orthogonal projection 
on the plane of T is that lozenge. 

If the tiling has no holes, the node of C we are finding ourselves at is independent of the way we traced 
around the tiles to get there, and what results is a lattice surface in C whose lattice square faces are in one 
to one correspondence with the tiles. 

When a hole of non-zero charge is present this ceases to be true. To illustrate this, suppose we have 
a left-pointing triangular hole of side two in our tiling (see Figure 5.1). Then as we trace its boundary 
counterclockwise, the traversal of each lattice segment agrees with its orientation. Thus, each complete 
turn results in six upward steps along lattice segments of £, and leaves us at the node of C which is two 
cube-body-diagonals higher than where we started (this is illustrated, from three different viewpoints, in 
Figure 5.2). For a similar right-pointing hole we would end up two such diagonals lower. 

In general, if we loop once counterclockwise along a closed walk through the lozenges of a tiling, the 
ending point of the corresponding walk on C is on the same vertical as its starting point, but a distance of 
below it, where q is the total charge of the holes we looped around (by a negative distance below we 

mean the absolute value of that distance above). This makes it impossible to get a single-sheeted lifting 
surface with no extra boundary in addition to the the lifting of the boundary of the region with holes. 

Nevertheless, we can construct a multi-sheeted surface with no additional boundary that lifts any tiling 
with holes, as follows. Suppose we have a lattice region TZonT with a finite number of holes, and let T be 
a lozenge tiling of it. For each hole, consider a fixed lattice path cut (independent of the tiling T) from it 
to the boundary of the region TZ] let V be the family of lattice paths formed by them. Whenever a step of 
a lattice path of V crosses a lozenge of T, remove that lozenge from TZ; let TZo be the region obtained from 
TZ by removing the union Tq of all such lozenges of T. Now regard T\Tq as being a tiling of TZq in which 
V n TZq is part of the boundary. Then T\Tq lifts to a single-sheeted lattice surface^ St in £. The union^ 

is a multi-sheeted surface having holes above all the lozenges in 7?. \ 7?.o. Define St to be the multi- 
sheeted surface obtained from U by filling in these holes with the missing square faces of the lattice £ 

^ We assume here that the paths in V are disjoint. Unless the holes are packed very close (which, since we are interested 
in scaling limits, will not be the case for us), V can be chosen so. 

® For a surface S (single- or multi-sheeted) embedded in and for c e M we denote by 5 -|- c the translation of S by the 
vector (0, 0, c). 
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Figure 5.1. A tiling with a triangular hole of size two. 






Figure 5.2. Three views of St for the tiling T above. 



above the lozenges in 7?. \ TZq. One readily sees that St is independent of the family of cuts V. The 
surface corresponding to the tiling in Figure 5.1 is illustrated in Figure 5.3. An instance with three holes 
is pictured in Figures 5.4 and 5.5. 

The detailed definition of the joint correlation a) is the following (sec [3]). For j — 1, . . . ,k, let Qj be 
either a lozenge-hole or a lattice triangular hole of side two. 

It is enough to define ^(Qi, . . . , Q^) when q = Y^'^=i'^^iQj) — (^^'^ other case reduces to this by 
reflection across a vertical lattice line). Our definition is inductive on q: 

(i). If q — 0, let N be large enough so that the lattice rhombus of side N centered at the origin encloses 
all Qj's, and denote by Tn the torus obtained from this large lattice rhombus by identifying its opposite 
sides. Set^ 

w(Qi, • ■ • ,Qfc) := lim . (5.1) 



N- 



{a). If g > 0, define 



u}{Qi, ...,Qk 



hm i?9(2;(Qi,...,Qfe,IF«,o) 



(5.2) 



The above limits exist by Proposition 2.1. 

Let Ai, . . . , Afc be fixed lattice triangular holes of side two, and let i be a fixed lozenge position. Assume 
q = X]j=i ch(Aj) = 0. In the limit measure of the uniform measures on the tori T/y \ Ai U . . . U A^., the 
probability that L is occupied by a lozenge in a random tiling is 



Prob {L is occupied} = lim 



M(TA,\LUAiU---UAfc) _ a;(L,Ai,...,Afe) 



N^oc M(r/v \ Ai U • • • U Afc) 



w(Ai 



,Afc) 



(5.3) 



^ M(7?.) denotes the number of lozenge tilings of the lattice region TZ. 
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Figure 5.3. (a) One of the two connected components of St for the tihng in Figure 5.1. (b) The fuU St- 



This expression for the probabihty of L being occupied holds in fact for general q. Indeed, suppose this 
has been established for total charges < q. Then the probability that L is occupied in a random tiling with 
the extra hole Wji q in addition to Ai, . . . , A^; is 

a)(X,Ai,...,Afc,Wfl,o) 

u{Ai,...,Ak,WR,o) 

However, in the limit R ^ oo this is, by (5.2), the same as the fraction on the right hand side of (5.3), and 
our statement is proved by induction. 

Suppose TZ is a bounded simply connected lattice region on T. Then each tiling of TZ lifts to a single- 
sheeted lattice surface on C, and we can define the average lifting surface S'av by taking the arithmetic 
mean of the finitely many heights of these surfaces above each node of T. As pointed out in [CLP] , if u 
and V are nearest neighbors in T so that the lattice segment between them is oriented from u to v (see the 
second paragraph of this section) , then 

SM - SM = 7f (1 ~ 3p(L™)) , (5.4) 

where L^v is the lozenge location whose short diagonal is uv (the factor multiplying the parenthesis on the 
right hand side arises because the traversal of each unit segment in T results in a change of height on the 
lifting surface of one third of a body diagonal of a lattice cube of £). Thus, the height of at any node 
of T can be obtained by taking cumulative sums of lozenge occupation probabilities. 

The regions we are concerned with — complements of finite unions of disjoint lattice triangular holes of 
side 2 — have an infinite set of lozenge tilings, so we cannot use the arithmetic mean as the definition of 
the average of the surfaces their tilings lift to. However, since we know the lozenge placement probabilities 
are given by (5.3), we can turn (5.4) around and use it to define this average surface. 

More precisely, consider for each triangular hole a lattice path in T from a point on its boundary to 
infinity. Let V be the union of these lattice paths; assume they are disjoint. Define to be the lattice 
surface on C satisfying 

Sav{v) - S'av(u) ^ 3p{Luv)) , (5.5) 

for any two nearest neighbors u and v of T for which the segment between them is oriented from u to v 
and is not crossed by any path of V. Define the average lifting surface of the tilings of the complement of 
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Figure 5.4. (a) A tiling T with three holes, (b) A view of St- 

the holes by 

Saw = S^-v + 71—=. 

nez V 2 

This definition is readily seen to be independent of the family of cuts V . 

We note that in the case when there are no holes, the average surface can be defined in terms of 
the translation invariant ergodic measure on the hexagonal lattice (the lattice whose dimer coverings are 
equivalent to lozenge tilings of the triangular lattice), which follows by a general result of Sheffield (see 
[14] [13]) to be unique and given by a limit of uniform measures on tori. Due to the presence of holes our 
setting does not seem to fit that context. 

The helicoid -ff (a, 6; c) is the surface whose parametric equations in Cartesian coordinates are 

X = a -\- p cos 6 

y = b + psmd, — oo < p, 6' < (X). (5-6) 

z^ce 

The half helicoid H~^{a, b; c) is obtained by restricting the range of p in the above to (0, oo). The dotted 
helicoid H{a, b; c) is H{a, b; c) minus the vertical axis x = a,y = b. 
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Figure 5.5. Two views of St for the tiling T in Figure 5.4(a). 
For positive integers s define the s-refined half hclicoid by 

H+ia, b; c) 1) H+{a, b;c) + ^. (5.7) 

Define the s-refined dotted helicoid by 



Hs{a,b;c):^\\H{a,b;c) + ^. (5.8) 
^ s 



Note that for s G Z, the fibers of H^{a,b; sc) above each point u in the xy coordinate plane are of the 
form f{u) + 27rcZ. Thus, given Hf\ai, bf, Sic), s; G Z, i = 1, . . . , fc, one can define their sum 

S := i7+(ai,5i;sic) H \- H^^{ak,bk; Skc) 

by defining the fiber of S above u to be 

hiu) + --- + fk{u) + 2TrcZ. (5.9) 
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Figure 5.6. (a) H+{0, 0; 2) + H+{0, 1; -1). (b) H+{0, 0; 2) + H+{0, 1; -1) is 
obtained by superimposing two copies of the surface (a). 

We can define the sum of dotted hehcoids analogously, using that the fibers of Hs{a, b; sc) above u are of 
the form f{u) + ttcZ. Figure 5.6(a) illustrates -ff^(0, 0; 2) + H^{1, 0; —1). Note that a point on the surface 
making a complete counterclockwise turn around the "spiral stairwell" on the left ends up two levels higher; 
a similar turning around the other spiral takes the point one level lower. Figure 5.6(b) illustrates the effect 
of doubling the refinement indices in a sum of refined half helicoids: the new surface is the union of the 
original surface and a suitable vertical translate of it. A depiction of a sum of three helicoids is given in 
Figure 5.7. 

Let Tr be the lattice obtained from the triangular lattice T by a homothety around the origin of factor . 

Theorem 5.1. Let x[^^ , . . . , Xm'^ , y[^'^ , . . . , ym'^ , z[^'' , ■ • ■ , ^n^'' and w[^'^ , ■ ■ ■ , wif'^ be sequences of integers 
so that limfl,^oo /R = Xi, limi^^oo v\^^ / R = Vr, limfl-^oo z^^^ / R = Zj and limfl^oo w^^'' /R = Wj for 
1 < i < m and 1 < J < n. Assume the {xi, yi) 's and (zj, Wj) 's are all distinct. 

Let -Elj, . . . , LJ^^^ and , . . . , , 3\l~q, be multiholes on Tr, and let the lists and hj have lengths 
Si and tj, respectively. Let 

be the average lifting surface of the tilings of the complement of these multiholes on Tr . 
Then, as R oo, RS^^ converges to the sum of refined helicoids 



i=l 



Si \ Xi^yi, - 



3st 
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i=i 



3^, 
V2tt 

it 



V27I 



(5.10) 



Figure 5.7. (a) The sum of the refined half hehcoids H+{-l, 0; -2), H+{0, 0; 3) and H+{1, 0; 4). (b) Two 
copies of the latter give 0; -2) + i76^(0, 0; 3) + iJ8^(l, 0; 4). 



For any bounded set B and any open set U containing (xi, j/i), . . . , {xm, Vm) cmd (zi, wi), . . . , (z„, Wn), the 
convergence is uniform on B\U. 

Proof. In addition to our 60° oblique coordinate system O in the plane of the unit triangular lattice T, 
consider also a Cartesian system of coordinates C having the origin at the node of T just below the origin 
of O , the X-axis in the polar direction and the y-axis in the polar direction 7r/2. One readily sees that 

the nodes of T have C-coordinates a + b even, and that the midpoint of the segment 

connecting the nodes {a^,b^^ and (a^, {b + 2)i^ of T has O-coordinates ^). 

We will use Cartesian coordinates to specify points at which S"^ is evaluated, and oblique coordinates 
for the lozenges whose placement probabilities come up. 

26 



We have by (5.5) 



' i? 2 ' i? 2 i I i? 2 ' ii 2 i V2 i 2ii ' 2R 



where pi is given by (2.12). The second equahty holds because placement probabilities are clearly invari- 
ant under scaling. Analogous equations relate the change of RS^^ when moving a distance of ^ from 

in the polar directions § + fcf , fc e Z. 



The asymptotics of pi [x\^ ,yQ 'j — p2 \yQ j follows by (3.1) and (1.2). One obtains the asymp- 
totics of pi {xq^'' , Uo^^^ —ps (^Xq^^ , Uo^^^ similarly from (1.3). Adding up the two and using pi +p2 +P3 = 1 



yields 

1 _ (r^R) - / V s- x^~x,+y^~ 

_ y^^ xq - Zj +yo- Wj I J: ^ ^ i_ 

^ ^ (xo - Zj)^ + {xo - Zj){yo - Wj) + {yo - Wjf j R ° \Rj ' ^.5^2) 

with the implicit constant uniform for {xo,yo) £ B\U. Similar expressions follow for the asymptotics of 
P2 and P3. 

This allows us to deduce that lim^^oo RSj^ exists. Indeed, the height of RS^^ at any node of 
can be obtained by taking cumulative sums of (5.11) and its two analogs. By (5.12) and its analogs, the 
resulting error terms have a zero effect in the limit R ^ 00: even though they do add up, and there are 
more and more of them as R gets large, since the individual errors are small and the factor is present, 
the combined contribution of the error terms is small. 

Furthermore, we can find the gradient of the limiting surface <S by dividing the left hand side of (5.11) by 
the distance between the two points where the heights of RSj^ are examined — and letting ii — > 00. 
This yields 

dS( VS l\ L/v^ "^{XQ-Xj) " "^{XQ-Zj) 

dy 2 '^'2) ^^'^L.'^^-^a^o-x.r + liyo-mr U^o - zj)^ + \{yo - wj) 

for (xo, yo) different from (xi, yi), . . . , {xm, ym) and {zi,wi), . . . , (z„, Wn)- After a change of variables this 
becomes 

dS I ^ J (XQ - Xj) -A (Xq - Zj) I /rioN 

^y(X0,2/0)- jZ.^'(^o_^.)2 + (y^_y.)2 Z.*^(^o_^.)2 + (y^_^.)2 j- 15-1"^) 

The analog of (5.11) corresponding to moving one unit in the polar direction ^ yields by a similar calcu- 
lation the directional derivative of <S along this polar direction. Writing djdx as the appropriate linear 
combination of this directional derivative and d/dy arrive at 

dS I ^ _ J (i/o - Vi) v^, (?/o - Wj) I 

dx - j A. (^^ _ ^.)2 + _ y.y - ^,)2 + (yo - Wj)^ J • ^'-'^^ 
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On the other hand, the gradient at (x, y) of the half hehcoid H~^{a, b; c) — and hence of any refinement 
of it — is (- (j,_a)2^+fa_b)2 , (^_a)2+(y-b)2 ^ ■ Thus, if H is the sum of refined heUcoids 



H 



i — 1 _7 — 1 



we have that the gradients of H and S agree on \ {(xi, yi), . . . , {zn, Wn)}- 

By construction, each lifting surface of a tiling of the complement of our multiholes has fibers of the form 
c + -^Z- Thus so does RSj^ for all R, and also the limit surface S. The same is true for each summand 
in (5.15), and therefore for H. Then the surface Sq = S — H is well defined by (5.9). By construction, as 
u approaches {xi,yi) the fiber of <S above u approaches + It follows from the definition of H that 
the same is true for the fibers of H. This, together with the fact that its gradient is zero, implies that <So 
is the surface + So S = H, justifying the first expression in (5.10). The second expression follows 

using the elementary fact that H^{a, b; c) U (i?+(a, h; fi) + f ) = H{a, b; c) (an instance of this can be seen 
in Figure 5.6(a)). The stated uniformity of the convergence is implied by the uniformity of the error terms 
in (5.12) and its analogs. □ 



6. Physical interpretation 



The results of [2] and [3] already show a close connection between random lozenge tilings with holes and 
electrostatics: the joint correlation of holes is given, in the scaling limit, by the electrostatic energy of the 
corresponding system of electrical charges. Theorem 1.1 strengthens this connection by showing that the 
electrostatic field itself can be viewed as the scaling limit of the discrete field of the average orientation of 
lozenges in a random tiling of the complement of the holes. 

Following [2, §2], suppose one pictures the two dimensional universe as being spatially quantized con- 
sisting of a very fine lattice of triangular quanta of space. 

Associate the complement of the holes with the vacuiim of empty space. The quantum fluctuations of the 
vacuum cause virtual electron-virtual positron pairs to be cciaselessly created and annihilated. Associate 
virtual electrons with left-pointing unit triangles and virtual positrons with right-pointing unit triangles. 
Then a lozenge corresponds to the annihilation of a virtual electron- virtual positron pair, and a lozenge 
tiling to one way for all virtual electron-virtual positron pairs to annihilate. 

Averaging over all possible tilings arises then naturally as corresponding to the Feynman sum over all 
possible ways the virtual electrons and virtual positrons can annihilate in pairs. Our result states that 
if the distances between the hole-charges and the distances between them and the point where we are 
measuring the F-field are much larger than the lattice spacing of the quantas of space, then this field is 
almost exactly the same as the electric field. Prom this perspective, one sees as the result of an exact 
calculation how the (two dimensional) electric field emerges from the quantum fluctuations of the vacuum. 
This view provides an electric field which very closely approximates the Coulomb field, but is not exactly 
equal to it: by Theorem 1.1, the latter occurs only in the limit when the size of the quantas of space 
approaches zero. The discrepancy gets noticeable only at very small distances. This is in agreement with 
the expectation of some physicists that linear superposition for the electric field might break down in the 
subatomic domain (see e.g. Jackson [11, §1.3]). 

In quantum mechanics, when a virtual electron and a virtual positron annihilate, two photons result, 
which are sent away in opposite directions. The unit contribution to F(e) of each lozenge covering e in 
a tiling could be regarded as corresponding to one of these two photons (this assumes that the common 
direction of the two ejected photons is, on average, the straight line connecting the annihilating virtual 
particles). Then Theorem 1.1 states that in the scaling limit, the thus generated "jet" of photons at e 
averages precisely to the electrostatic field. This could then be viewed as a "mechanism" for the coming 
about of the electrostatic field. 

The discussion in the previous paragraph leaves the second photon unaccounted for. However, underlying 
in our considerations there is a second discrete vector field we have not mentioned yet: the one given by 
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the average orientation of the lozenge that covers any fixed n(7/it-pointing unit triangle! Denote it by F'. 
An argument similar to the one that proved Theorem 1.1 shows that F' = — F in the scaling limit. 

We conclude by mentioning another, quite different way to define a vector field via random tilings with 
holes. Let Qi, . . . , Qn be a fixed collection of holes on the unit triangular lattice. For any x,y,a, (3 € Z 



define 



J. / N _ ^ f ^ iEx+a,y+fj, Ql, - ■ ■ , Qn) _ \ 



(i.e., the hole E^^y plays the role of a "test charge," and the effects of its displacement are recorded). Then 
it can be shown (details will appear in a separate paper) that there exists a vector field T so that in the 
scaling limit Ta^f}{x, y) is the orthogonal projection of T(x, y) onto the vector (a, [3). Furthermore, up to a 
constant multiple, the field T turns out to be the same as F. Thus, unlike in physics where the electric field 
is defined by means of a test charge, in our model there are two different ways to define the corresponding 
field, and one of them does not use test charges. 

We note that the agreement in the previous paragraph is not a matter of course, as it does not hold for 
instance on the critical Fisher lattice — the planar lattice of equilateral triangles and regular dodecagons 
where the edges of the triangles have weight 1, and the inter-triangle edges have weight a/S (this follows 
from not yet published joint work with David Wilson). 
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